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We study the low-energy dynamics of an acoustic black hole near the sonic 
horizon. For the experimental test of black hole evaporation in the laboratory, 
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5— i . the decay rate (greybody factor) of the acoustic black hole (sonic hole) can be 



calculated by the usual low-energy perturbation method. As a consequence, 

we obtain the decay rate of the sonic horizon from the absorption and the 

reflection coefficients. Moreover, we show that the thermal emission from the 

sonic horizon is only proportional to a control parameter which describes the 

velocity of the fluid. 
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One of the most attractive aspects of general relativity is that it predicts the existence 
of black holes. Classically, they can absorb all sort of matter through the extremely strong 
gravitational field, and even light cannot escape from them. For this reason, they are black 
to all observers outside them. Their quantum mechanical properties are quite different from 
their classical ones. The surprising prediction that quantum black holes are no longer black 
was announced in 1974 [1], thereafter the quantum mechanical properties of them remain 
fascinating but problematic to theoretical and experimental physicists. From the viewpoint 
of quantum field theory in a curved spacetime background, they can emit thermal radiation 
proportional to the surface gravity at their horizons. This theoretical result, called "Hawking 
radiation" , used to be asked if it is possible to verify it from experimental study within the 
framework of astrophysical phenomenology. 

In 1981 Unruh proposed a possibility of an experimental test for the black hole evap- 
oration [2], which opened a new branch of black hole physics [3-12]. He considered an 
irrotational inviscid fluid and sound wave propagating in the medium. According to Ref. 
[2], the linearized small perturbation of wave equations in the fluid leads to an equation 
of motion for a massless scalar field in the geometrical background that is similar to the 
black hole metric. There exists a trapped region that the sound wave cannot escape as the 
fluid velocity exceeds the sound velocity (hypersonic flow), which is called a "sonic horizon" 
when the fluid velocity is equals to that of the sound wave. This hypersonic flow system 
reveals more aspects of the acoustic black hole with the help of the analogy of black hole 
thermodynamics. By quantizing the sound field (phonons), the sonic hole (dumb hole) can 
emit sound waves with the thermal temperature given by 

nn & 

T = — 

2vr 

where h is a Planck constant and v r is a fluid velocity along the radial direction. 

In this paper, we shall investigate the decay rate of sonic black holes from the near- 
horizon low-energy dynamics, which is expected to give an experimental suggestion of de- 
tecting the thermal radiation. The low-energy perturbation analysis of phonons leads to an 
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exact derivation of the absorption cross section and decay rate of the sonic horizon. As a 
consequence, those results depend only upon a control parameter that describes the veloc- 
ity of the given fluid associated with the surface gravity on the horizon, which agrees with 
Unruh's prediction of thermal radiation. 

Let us briefly review the hypersonic flow and acoustic black holes following Ref. [2] . The 
irrotational inviscid fluid flow is described by the following three equations, 



Vxt?=0, 
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dt 
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p + V • (pv) = 0, 
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where v is a velocity of the flow, p is a pressure that is a function of the density p, and 
<3> is an external potential. Note that the first equation describes an irrotational flow, the 
second and the third ones are Eulerian and continuity equations, respectively. These can 
be linearized in the vicinity of some mean flow solution with £ = £o + £ and ip = ipo + ip 
redefining the fields as VG(£) = Vp/p, p = e^, and v = V0, where the barred parameters 
stand for the small fluctuations around the vacuum solutions with a zero label. Combining 
these linearized equations gives the massless scalar field equation describing the sound wave 
in the curved spacetime background 
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with the background metric, 
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where i,j = 1,2,3 and c is a local velocity of sound defined as c 2 = G'(£). Following 
Ref. [2], we assume c to be constant for convenience. 1 . Provided that the metric, (4), 



This is assumed to be constant since the role of sound in the fluid is similar to that of light in 
gravity, aside from the easy treatment of the remaining calculations. 
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is spherically symmetric and stationary, we can change the time coordinate t into r = 
t + f dr{v^{r) I '(c 2 — t>Q 2 (r)), which yields a diagonalized metric, 
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The coordinate transformation does not affect the remaining physical result although the 
metric is a coordinate-dependent quantity. We, therefore, can make it deform for conve- 
nience. Of course, if the metric had been transformed in the scalar field equation of Eq. 
(3), then the same diagonalized metric would have been obtained. The metric (5) has a 
sonic horizon at c 2 = v^ 2 , which coincides with the sonic fluid flow while the hypersonic flow 
(c 2 < Vq 2 ) prevents the sound wave from propagating and the sound wave eventually does 
not escape from the horizon. 

If the flow smoothly exceeds the velocity of sound at the sonic horizon r = r sh , the 
velocity can be expanded as Vq(t) = — c + a(r — r s h) + 0(r — r s h) 2 , which yields (up to the 
first order of r) 



c 2 - v r G 2 ps lacir - r sh ), 



(6) 



where a is a control parameter associated with the velocity of fluid defined as (V • v)\ r=rsh . 
Note that the parameter a is related to the surface gravity up to some proportional factor, 
which will be shown later. The resulting metric, therefore, becomes 

cdr 2 
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(ds) =- 



-2ac(r - r sh )dr + 
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which is similar to the near-horizon behavior of the Schwarzschild black hole as discussed in 
Ref. [2]. 

Our starting point is the field equation of the sound wave, Eq. (3), under the near- 
horizon metric background of Eq. (7). Using the separation of variables, ijj(T,r,9,(p) = 
u{r-)Q(6)e i{mLp - wr \ the radial equation of Eq. (3) is written as 



(r - r sh )d 2 r u{r) + 



1 + -(r - r sh ) 
r 



d r u(r) 
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u(r) = 0. (8) 



What we wish to do is to approximate the wave equation in the Schwarzschild geometry 
by another unphysical wave equation that has the same behavior near the sonic horizon. 
Even though the unphysical wave equation has the tremendous advantage that it is an- 
alytically solvable, we will work within the first order of the approximation, because the 
near-horizon limit of the Schwarzschild metric gives the simplified Eqs. (7) and (8). 
From the change of variables, < z — 1 — r s h/r < 1, the radial equation (8) becomes 
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z(l — z)d z u(z) + d z u(z) 
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u(z) = 0, 



(9) 



z(l-z) 

where A = c£(£ + l)/2ar s h and B = u 2 /4a 2 . Now let us define u(z) = z a (l — z) f5 h(z) in 
order to remove singularities at z = and z = 1, which leads to 

z(l - z)d 2 z h(z) + [1 + 2a - 2(a + (3)z\d z h(z) - [A + (a + (3)(a + (3 - l)]h(z) = 0, (10) 

and the parameters a and (3 are determined by a = ±iy^B and (3 = l±Vl — B, respectively. 
The solution of Eq. (10) is simply given by the hypergeometric function F as 

u(z) = C lZ a (l - zfF ^~ + a + (3- ivi-4A, ~ + a + (3 + \^f\ - AA; 1 + 2a; z^j 
+ C 2 z~ a (l - zfF - a + (3 - ^Vl - 4A, ~- a + /3+ - AA; 1 - 2a; z^j . (11) 

Note that since this solution is symmetric for changing a into —a, we take the plus sign of 
a for convenience. In z — > limit, Eq. (11) becomes 
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where C out = C\ and C in = C 2 . At this stage, the boundary condition can be imposed 
with an appropriate physical situation, and here we take C out = 0. The configuration of the 
boundary condition that we are considering is shown in the FIG. 1. Using the z — > 1 — z 
transformation of hypergeometric function [13], 

x r(n)F(n - I - m) 
F(l, m; n; z) = _ _ F(l, m; I + m - n + 1; 1 - z) 

+ (1 - z )n-;-m £MSj + m ~ n ) ^( n _ I. n _ m . n _ I _ m + 1; 1 _ z \ (13) 

T(l)T(m) 



sonic horizon 




FIG. 1. Configuration of the acoustic black holes and the boundary condition 



the far- region solution of Eq. (11) for z — > 1 is given by 

Ufar(r) ~ Q 



r(l - 2a)r(2 - 2/3)/ 



sh -p 



+ 
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r(l-2a)r(2/?-2)r,y ,_ 2 



(14) 



r(-±-a + P- I VI - 4 A) T(-\-a + {3+ - A A) 
This far region solution is symmetric for changing j3 into 2-/3, which means that we can 
take the plus sign of (3 for convenience. Since (3 = 1 + y/Aa 2 - uo 2 /2a, we consider the real 
value case of (3, which signifies that the control parameter of the fluid velocity, a, is much 
greater than the energy of sound wave, 4a 2 > u 2 . This is nothing but a low-energy situation 
that we are considering. 

On the other hand, the asymptotic behavior of Eq. (8) produces the wave equation at 
boundary, 
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r 2 <9> B (r) + 3rd r u B (r) + ^b(0 = 0, (15) 
which can be solved by 

u B (r) = A ir ~P + A 2 r-P- 2 . (16) 



The coefficients, C m and Ai/A^), can be exactly matched by comparing Eq. (14) with 
Eq. (16). Note that f3 — 2 at low frequency region since 4a 2 >> oo 2 . The "ingoing" and 
"outgoing" waves (denoted by A in and A out , respectively) of Eq. (16) can be decomposed 
by introducing the redefinition of amplitudes, 



Ax = iL(A in - A out ), 
A 2 = A- m + A out , 



(17) 



which leads to 
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out 
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(18) 



where L is a positive dimensionless constant, which will be taken to be independent of 
the energy uj [14,16]. Note that L can be chosen so that the absorption cross-section has 
the value of the area of the black hole as u — > [14,17] or it has the usual value of the 
Hawking temperature [16]. However, in this approach, we shall take into account a different 
requirement of determining a positive dimensionless constant L. The fact that A + 1Z — 1 
can determine the dimensionless constant L in our configuration shown in FIG. 1. 

The absorption (A) and the reflection (71) coefficients can be defined by the ratio of 
"ingoing" and "outgoing" fluxes, 
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respectively. It is easy to derive the flux from Eq. (8) and it becomes 

2-7T 

T = —r 2 g rr (u*(r)d r u(r) — u(r)d r u*(r)). 
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(20) 



The reflection coefficient can be expressed in terms of the matched coefficients in Eqs. (14) 
and (16) as 
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For I = and uo — > 0, the \A 2 \ 2 term in Eq. (21) is evaluated by 
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' "' ' mi T(cu a -a)T(l + uj a -a)T(cu a + a)T(l + u a + a) 

~ |Cin| r sh , (22) 

by using T(z)T(l — z) = n/ sinnz. This approximation means that the Compton wavelength 
of the scattered wave is much greater than the sonic black hole size. Furthermore, the 
dominant contribution of s-mode is the case of vanishing angular potential^ = 0) of the 
whole propagation modes [18]. Similarly, the \Ai\ 2 term can be calculated as \Ai\ 2 m \C- m \ 2 rj h 
in the same limit used in Eq. (22). Provided that this is compared with Eq. (22), \Ai\ 2 is 
negligible for the small sonic black hole horizon. 

On the other hand, the non-trivial contribution of the reflection coefficient arises from 
the term of A\A X - A\A 2 . Using r(z)r(l - z) — ir/ sin tcz, we obtain 

A* 2 A t - A\A 2 w l r ^ UJ \ C ^\ (23) 

2a 

From Eqs. (22) and (23), the reflection coefficient is simply expressed by 

TZ e -° « , 24 

which is expected to have a range with < 1Z < 1 by definition. This fact restricts the 
frequency of u to < uj < 2aL. The absorption coefficient can be obtained by the reflection 
coefficient since we have 1 — TZ = A, which produces A e=0 = 2uo/(2aL + uj). From the 
definition of absorption coefficient in Eq. (19), we also obtain A i=0 = uj/2L 2 (2aL + a;). 
Comparing these both relations determines the positive dimensionless constant, L = 1/2. 
As a consequence, the reflection and the absorption coefficients are expanded as 

a a 

and their behaviors are plotted in the FIG. 2. Note that the critical value of a c = 2uj exists, 
which gives "the total absorption" (A i=0 = 1) and "the no reflection" (1Z i=0 = 0). This 
signifies that we have to adjust the relation between the velocity and the energy of phonons, 
when we control the velocity of the fluid in the experimental system. 
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FIG. 2. Plots of the reflection and the absorption coefficients for lo = 1. The valid range of a 
is a > 1 for the fixed u = 1. 



On the other hand, a relation between the vacuum expectation value of the number 
operator and the reflection coefficient can be given in the form [15,16], 
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Finally, the Hawking temperature is defined by 



< 0|iV|0 >= 



exp(u/T H )-V 

and combining this with Eq. (25) leads to the Hawking temperature as 



(27) 



H\oj->0 



a 
2' 



(28) 



This result coincides with the well-known result of the Hawking temperature in that it 
depends only upon the flow divergence (a) passing through the sonic horizon surface at a 
certain time, which is closely related to surface gravity on the horizon. The absorption cross 
section for the dominant s-wave mode is given by = A/uj = 2/(a + u>) with the help of 
Eq. (25), and the decay rate of the sonic horizon can be easily calculated by 



T^=0 

decay 



a 



abs 



e ^/T H _ I 
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The system we are considering is only for the low-energy frequency of the phonon field. 
However, at high frequencies, there exist a kind of backscattering from the geometry of 
the acoustic black hole. More precisely, this is equivalent to the case of uj 2 > 4a 2 in our 
approach, which yields an imaginary value of (5. A straightforward calculation following a 
similar procedure at high frequencies produces the exact results for the reflection and the 
absorption coefficients, 

COsh -KUJa + COS jlTC 



COsh 71LU+ + COS jfTT ' 

:osh 7tc<j+ — cosh iru>' c 

COsh 7TCJ+ + COS jglV 



cosh 7Tc<j+ — cosh 1X0J n , , 

**»> = ^ ' + I Z^T > (30) 



respectively. Note that to^ — u/a ± \Juj 2 — Aa? /a and j e = yl — ^^-^(^ + 1) and it is easy 
to show that TZuj» + Aj>> = 1- Since oj+ ~ 2w/a and cj^ w at high frequency region, 
the reflection (71) and the absorption coefficients (A) are 

1 + cos jeii 



cosh 27ru;a + cosj^n 
cosh 27ru;a — 1 

"4a;>> 



a;>> 



0, 

1. (31) 



cosh 27ra;a + cosj^n 

Note that this behavior is nothing but a scattering problem in the flat spacetime background. 
In the tortoise coordinate system, r* = f dr/a(r), one may read off the effective potential 
given by V e s(r) = u 2 — Ve(r) from the Wheeler- Regge equation, where 



V e (r) = a(r) 



d r a(r) | 1(1+1) 



(32) 



and cr(r) is a metric function in the metric, (ds) 2 = —a(r)dt 2 + a~ l (r-)dr 2 + d 2 Q. It is easily 
shown that the effective potential of the wave equation becomes trivial as V e s(r) = uj 2 at 
high frequencies. Therefore, the Wheeler-Regge equation describes the free field equation 
with the constant effective potential of flat spacetime. Furthermore, the decay rate at high 
frequencies from Eq. (29) vanishes as u goes to infinity, which is compatible with the 
statements at high frequencies as discussed before. 



10 



ACKNOWLEDGMENTS 



J.J.O. would like to thank Edwin J. Son, Jaehoon Kim, and Hyunjoo Lee for exciting 
and valuable discussions and he is grateful to Yoonbai Kim, Hyeonjoon Shin, and Gungwon 
Kang for a warm hospitality while he was staying at SKKU and KIAS as a visiting scholar. 
He also wishes to thank Russell J. Smith for helpful comments on this paper. This work 
was supported by the Basic Science Research Program of the Korean Science & Engineering 
Foundation Grant. At the final stage of this work, J.J.O. was supported by the Post-doctoral 
Fellowship Program of Korea Science & Engineering Foundation (KOSEF). 



11 



REFERENCES 

[1] S. W. Hawking, Nature, 248 (1974) 30. 

[2] W. G. Unruh, Phys. Rev. Lett. 46 (1981) 1351. 

[3] M. Novello, M. Visser, and G. Volovik, Artificial black holes, World Scientific, River 
Edge, USA (2002). 

[4] S. Corley and T. Jacobson, Phys. Rev. D57 (1998) 6269. 

[5] T. Jacobson and G. E. Volovik, Phys. Rev. D58 (1998) 064021. 

[6] T. Jacobson and D. Mattingly, Phys. Rev. D61 (1999) 024017. 

[7] W. G. Unruh, Phys. Rev. D51 (1995) 2827. 

[8] W. G. Unruh and R. Schwtzhold, Phys. Rev. D68 (2003) 024008. 

[9] M. Visser, Class. Quantum Grav. 15 (1998) 1767. 

[10] S. Liberati, S. Sonego, and M. Visser, Class. Quantum Grav. 17 (2000) 2903. 

[11] C. Barceld, S. Liberati, and M. Visser, Class. Quantum Grav. 18 (2001) 1137. 

[12] R. Balbinot, S. Fagnocchi, and A. Fabbri, Quantum Effects in Acoustic Black Holes:The 
Backreaction, [gr-qc/0405098]; R. Balbinot, S. Fagnocchi, A. Fabbri, and G. P. Proco- 
pio, Backreaction in Acoustic Black Holes, [gr-qc/0405096]. 

[13] M. Abramowitz and I. A. Stegun, Handbook of Mathematical Functions (Dover Publi- 
cations Inc., New York, ninth printing, 1970). 

[14] D. Birmingham, I. Sachs, and S. Sen, Phys. Lett. B413 (1997) 281. 

[15] K. Ghoroku and A. L. Larsen, Phys. Lett. B328 (1994) 28. 

[16] W. T. Kim and J. J. Oh, Phys. Lett. B461 (1999) 189. 

[17] S. Das, G. Gibbons, and S. Mathur, Phys. Rev. Lett. 78 (1997) 417. 

12 



[18] W. T. Kim, J. J. Oh, and K. H. Yee, Phys. Rev. D66 (2002) 044017. 



13 



